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It  is  well  known  and  demonstrated  that  interference  spectroscopy  offers  capabilities  to  obtain  passive 
remote  optical  sensing  spectra  of  high  precision  and  also  achieves  economies  in  size,  cost,  and  ease  of 
deployment  compared  with  more  conventional  systems.  We  describe  the  development  of  a  near¬ 
ultraviolet  spatial  heterodyne  spectrometer  designed  for  remote  sensing  of  the  global  distribution  of  the 
hydroxyl  radical  OH  in  the  Earth’s  middle  atmosphere.  The  instrument,  known  as  SHIMMER  (Spatial 
Heterodyne  Imager  for  Mesospheric  Radicals),  is  expected  to  obtain  its  first  OH  measurement  from  space 
in  early  2002  from  the  Space  Shuttle.  ©  2002  Optical  Society  of  America 
OCIS  codes:  010.1290,  120.3180,  120.6200. 


1.  Introduction 

The  instrument  described  in  this  paper  draws  its 
primary  scientific  motivation  from  the  Middle  Atmo¬ 
sphere  High  Resolution  Spectrograph  Investigation 
(MAHRSI),  which  was  designed  specifically  to  make 
global  measurements  of  OH  in  the  middle  atmo¬ 
sphere  from  low  Earth  orbit.  The  importance  of  the 
hydroxyl  radical  OH  in  the  photochemistry  of  plane¬ 
tary  atmospheres  has  been  recognized  for  more  than 
50  years1;  nevertheless,  recent  satellite  observation 
with  high-resolution  ultraviolet  (UV)  spectroscopy 
has  resulted  in  a  number  of  unexpected  results.  The 
MAHRSI  measurements  provide  firm  evidence  that 
the  photochemistry  of  ozone  and  the  odd-hydrogen 
HO*  family  (H,  OH,  and  H02)  as  formulated  in  cur¬ 
rent  atmospheric  models  cannot  consistently  explain 
the  abundance  of  OH  throughout  the  middle 
atmosphere.2-4  In  addition,  comparisons  of  the 
MAHRSI  data  with  direct  measurements  of  meso¬ 
spheric  H20  by  the  Halogen  Occultation  Experiment 
(HALOE)  instrument  on  the  Upper  Atmosphere  Re- 
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search  Satellite  (UARS)  have  revealed  that  the  OH 
data  also  comprise  a  powerful  measure  of  the  distri¬ 
bution  of  mesospheric  water  vapor.2  A  layer  of  H20 
detected  near  an  altitude  of  65  km  violates  our  un¬ 
derstanding  of  the  processes  by  which  water  is  pro¬ 
duced  in  and  transported  through  the  middle 
atmosphere.  MAHRSI  data  also  led  to  the  discovery 
of  a  water-vapor  layer  near  the  Arctic  summer  me- 
sopause,  possibly  related  to  the  formation  of  high- 
altitude  (82-km)  polar  mesospheric  clouds;  it  has 
been  suggested5  that  the  increasing  occurrence  of 
such  clouds  is  a  vivid  indication  of  atmospheric 
change  caused  by  human  activity.  However,  to  our 
knowledge,  there  is  little  quantitative  knowledge  of 
why  polar  mesospheric  clouds  form  and  why  they 
vary. 

In  its  measurement  of  OH,  MAHRSI  observed  the 
UV  solar  resonant  fluorescent  band  of  OH  near 
308-nm  wavelength  with  a  resolving  power  of  approx¬ 
imately  15,000,®  which  was  high  enough  to  separate 
unambiguously  the  OH  emission  from  the  bright 
spectrally  complex  background  formed  by  sunlight 
Rayleigh  scattered  by  the  ambient  atmosphere.  The 
MAHRSI  results  and  the  atmospheric  science  ques¬ 
tions  they  raise  clearly  define  the  need  and  the  re¬ 
quirements  for  a  high-spectral-resolution  UV 
instrument  with  high  spatial  resolution  and  high 
throughput.  However,  the  weight  (138  kg),  size  ( 120 
cm  X  85  cm  X  23  cm),  complexity,  and  responsivity  of 
a  conventional  spectrograph  like  MAHRSI  make  it 
unsuitable  for  future  satellite  mission  opportunities. 

Using  an  innovative  optical  technique  called  spa¬ 
tial  heterodyne  spectroscopy  (SHS),7’8  we  developed 


1  March  2002  /  Vol.  41,  No.  7  /  APPLIED  OPTICS 


1343 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

MAR  2002  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2002  to  00-00-2002 

4.  TITLE  AND  SUBTITLE 

SHIMMER:  a  spatial  heterodyne  spectrometer  for  remote  sensing  of 
Earth’s  middle  atmosphere 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory ,4555  Overlook  Avenue 

SW, Washington, DC, 20375 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

10 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


an  instrument  called  the  Spatial  Heterodyne  Imager 
for  Mesospheric  Radicals,  or  SHIMMER,  which  ap¬ 
pears  capable  of  outperforming  MAHRSI  in  a  much 
smaller  and  lighter  package.  Extensive  laboratory 
testing  indicates  that  the  performance  of  SHIMMER 
is  near  theoretical  predictions.  Despite  its  small 
size  and  great  simplicity,  laboratory  calibrations 
demonstrate  that  the  responsivity  of  SHIMMER 
[(counts/s)/(kRayleigh/nm)]  is  more  than  10,000 
times  larger  than  MAHRSI  and  its  spectral  resolu¬ 
tion  (0.0058  nm)  is  3.5  times  better.  Its  mechanical 
volume  and  mass  are  approximately  17%  of  MAHR- 
SI’s.  A  proof-of-concept  test  flight  of  SHIMMER  is 
planned  from  the  middeck  of  the  Space  Shuttle  by  use 
of  the  high-quality  optical  window  in  the  orbiter  side 
hatch.  Environmental  tests  have  verified  that  the 
design  is  suitable  for  Shuttle  launch  environments. 

In  the  following  sections  we  describe  the  funda¬ 
mental  concepts,  design,  laboratory  evaluation,  his¬ 
tory,  and  current  status  of  SHIMMER. 

2.  Properties  of  Spatial  Heterodyne  Spectroscopy 
Instruments 

SHS  is  particularly  well  suited  to  UV  applications 
that  require  high  resolving  power  and  high  through¬ 
put  over  a  relatively  narrow  spectral  band.  In  gen¬ 
eral,  however,  a  spatial  heterodyne  spectrometer  can 
be  tuned  to  work  in  any  spectral  region  with  control 
over  the  trade-off  between  spectral  resolution  and 
bandwidth.  SHS  is  similar  to  Fourier-transform 
spectroscopy  (FTS),  with  three  primary  differences. 
First,  SHS  effectively  heterodynes  the  interferogram 
about  a  selected  reference  wavelength,  resulting  in 
easily  measurable  low-spatial-frequency  fringes  pro¬ 
duced  by  wavelengths  in  the  passband  of  the  instru¬ 
ment.  Second,  the  spatial  heterodyne  spectrometer 
has  no  moving  parts,  even  in  the  field-widened  mode, 
and  therefore  is  much  smaller,  lighter,  and  less  com¬ 
plex  mechanically  than  FTS  instruments  making 
similar  measurements.  Third,  SHS  uses  an  imag¬ 
ing  detector  to  record  a  fringe  pattern  localized  inside 
the  interferometer,  consequently  phase  errors  result¬ 
ing  from  optical  defects  can  be  easily  corrected  in  the 
data  reduction.  Compared  with  more  conventional 
grating  spectrometers  such  as  MAHRSI,  SHS 
achieves  orders  of  magnitude  higher  throughput  at 
similar  or  higher  spectral  resolution  in  a  much 
smaller  package.8  As  is  the  case  with  FTS  instru¬ 
ments,  zero,  one,  or  two  dimensions  of  spatial  infor¬ 
mation  can  be  obtained  by  appropriate  choice  of  the 
optics  prior  to  the  interferometer. 

A.  Basic  Spatial  Heterodyne  Spectroscopy  On-Axis 
Analysis 

The  basic  SHS  concept  is  depicted  in  the  Fig.  1,  where 
the  image  of  an  extended  scene  is  collimated  to  form 
a  beam  incident  on  the  beam  splitter.  Diffraction 
gratings  in  each  arm  return  the  light  to  the  beam 
splitter.  Lenses  L2  and  L3  in  the  exit  beam  image 
the  gratings  onto  a  position-sensitive  detector.  The 
gratings  are  positioned  so  that  they  are  the  same 
distance  from  the  beam  splitter  along  the  optical  axis 


Fig.  1.  Conceptual  diagram  of  a  field-widened  spatial  heterodyne 
spectrometer. 


and  are  set  at  an  angle  0L  such  that  rays  of  a  certain 
wavelength  incident  on  the  gratings  parallel  to  the 
optical  axis  retroreflect,  returning  parallel  to  the  op¬ 
tical  axis.  At  this  wavelength,  called  the  Littrow 
wavelength,  the  recombining  wave  fronts  exit  the  in¬ 
terferometer  completely  in  phase  with  each  other 
combining  constructively  at  all  points  across  the  ap¬ 
erture.  For  a  wavelength  slightly  different  from  the 
Littrow  wavelength,  the  diffracted  beams  return  at  a 
small  angle  to  the  optical  axis,  and  the  recombining 
wave  fronts  are  crossed  as  depicted  in  Fig.  1.  In  this 
case  the  phase  difference  between  the  two  wave 
fronts  changes  from  a  maximum  at  one  edge,  to  zero 
in  the  middle,  and  back  to  maximum  at  the  other 
edge,  resulting  in  a  Fizeau  fringe  pattern  whose  spa¬ 
tial  frequency  is  related  to  the  angle  at  which  the 
wave  fronts  cross. 

The  off-Littrow  angle  y,  at  which  incident  light 
parallel  to  the  optical  axis  of  wave  number  <r  exits  the 
interferometer,  is  given  by  the  grating  equation 

o-[sin  0L  +  sin(0L  -  y)]  =  m/d,  (1) 

where  m  is  the  order  of  diffraction,  0 L  is  the  Littrow 
angle,  and  1/d  is  the  grating  groove  density.  The 
spatial  frequency  of  Fizeau  fringes  is  related  to  the 
wave  number  of  light  by  the  equation 

fx  =  2a  sin  y  ~  4(<r  —  o-0)tan  0L,  (2) 

where  the  approximation  assumes  small  y,  and  ct0  is 
the  Littrow  wave  number  for  which  y  =  0.  For  axial 
rays,  the  intensity  at  the  fringe  localization  plane  for 
incident  spectral  density  B(a)  is  given  by 

I{x)  =  B(a)(  1  +  cos{2tt[4(ct  —  a0)x  tan(0L)]})do-, 

Jo 

(3) 
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where  x  is  measured  in  the  dispersion  plane  of  the 
gratings.  Note  that,  except  for  the  constant  term, 
I(x)  is  the  Fourier  transform  of  the  input  spectrum. 
The  fringe  localization  plane  is  at  the  gratings,  and 
therefore  the  gratings  are  imaged  directly  on  the  de¬ 
tector  to  produce  an  interferogram  that  is  then  Fou¬ 
rier  transformed  to  recover  the  input  spectrum.  The 
maximum  achievable  resolving  power  is  given  by 
R0  =  4  Wcr  sin  0L  where  W  is  the  width  of  the  gratings. 
It  can  be  shown  that  R0  is  simply  equal  to  the  total 
number  of  grating  grooves  imaged  onto  the  detector. 

The  nonaliased  passband  of  the  spatial  heterodyne 
spectrometer  is  limited  by  the  number  of  samples 
across  the  interferogram.  If  the  detector  has  N  pix¬ 
els  in  the  dispersion  plane  of  the  gratings  (the  plane 
of  Fig.  1),  N/ 2  nonaliased  spectral  resolution  ele¬ 
ments  can  be  measured,  resulting  in  a  passband  of 
Ncr/{2R).  To  achieve  this  maximum  bandwidth,  the 
optics  that  image  the  gratings  on  the  detector  must 
adequately  reproduce  N/2  fringes  at  the  detector. 

The  symmetry  of  the  cosine  interferogram  in  Eq. 
(3)  results  in  an  ambiguity  between  wave  numbers 
above  and  below  the  Littrow  wave  number  (a  ±  Act) 
as  their  fringe  frequencies  are  identical.  This  has 
the  effect  of  folding  the  spectrum  about  the  Littrow 
wave  number.  This  ambiguity  can  be  handled  in 
one  of  three  ways:  (1)  The  cosine  symmetry  can  be 
broken  by  the  introduction  of  a  small  tilt  to  one  of  the 
gratings  to  produce  low-spatial-frequency  fringes  in 
the  y  dimension  of  the  detector  in  addition  to  the 
wave-number-dependent  fringes  in  the  x  direction  il¬ 
lustrated  in  Eq.  (3).  With  this  y  tilt,  the  resulting 
fringe  patterns  are  rotated  clockwise  or  counterclock¬ 
wise  depending  on  whether  ct  is  above  or  below  the 
Littrow  wave  number.  A  two-dimensional  Fourier 
transform  distinguishes  between  these  fringe  pat¬ 
terns.  (2)  A  prefilter  can  be  used  to  block  light  on 
one  side  of  the  Littrow  wave  number.  (3)  The  folded 
spectrum  can  be  accounted  for  in  the  spectral  analy¬ 
sis.  The  pre filter  solution  to  this  problem  was  cho¬ 
sen  for  SHIMMER. 

B.  Off-Axis  Analysis 

When  off-axis  angles  are  included  in  the  grating 
equation,  it  can  be  shown7  that  the  solid-angle  field  of 
view  of  the  interferometer  without  field-widening 
prisms  is  characteristic  of  a  conventional  interferom¬ 
eter  (Fabry-Perot  or  FTS)  operating  at  the  same  re¬ 
solving  power  and  is  given  by  If  =  2rn/R0.  This  limit 
is  defined  as  the  field  of  view  at  which,  for  monochro¬ 
matic  light,  extreme  rays  produce  one  additional 
fringe  across  the  detector  compared  with  the  axial 
rays.  With  the  field  of  view  limited  to  Q,  a  mono¬ 
chromatic  interferogram  will  have  highest  fringe  con¬ 
trast  in  the  center  where  all  fringe  patterns  from  all 
angles  are  in  phase  and  a  slightly  decreasing  contrast 
toward  the  edges  as  the  fringes  become  increasingly 
out  of  phase.  Because  of  this  phase  difference  at  its 
edges,  the  envelope  of  the  interferogram  changes 
from  a  rectangular  function  to  a  rectangle  multiplied 
by  a  sine  function.7  The  resulting  instrumental  pro¬ 
file  is  broadened  slightly  and  resolution  is  reduced 


slightly  as  a  result.  A  further  increase  in  the  field  of 
view  will  result  in  reduced  fringe  contrast  at  the  ends 
of  the  interferogram  and  an  unacceptably  broadened 
instrument  function.  These  considerations  are  ex¬ 
actly  analogous  to  the  trade-off  between  field  of  view, 
fringe  amplitude,  contrast,  and  resolution  in  Michel- 
son  FTS.9 

The  field  of  view  can  be  greatly  increased  through 
use  of  the  field- widening  prisms  shown  in  Fig.  1. 
The  fixed  prisms  serve  to  rotate  the  images  of  the 
gratings  (from  a  geometrical  optics  point  of  view)  so 
that  they  appear  to  be  normal  to  the  optical  axis.  In 
this  configuration,  the  spatial-frequency  dependence 
with  an  off-axis  angle  is  minimized  for  monochro¬ 
matic  light,  and  fringe  patterns  with  nearly  constant 
contrast  across  the  entire  detector  are  produced  even 
when  the  field  of  view  at  the  gratings  is  increased  by 
an  order  of  magnitude  above  the  non-field-widened 
case.  Such  a  large  field  of  view  gives  SHS  an  enor¬ 
mous  throughput  advantage  over  grating  spectrom¬ 
eters  that  require  a  narrow  slit  to  achieve  high 
resolving  power. 

A  detailed  analysis  of  field  widening  shows  that  the 
field  of  view  of  the  interferometer  is  limited  by  either 
prism  spherical  aberration  for  systems  with  small 
Littrow  angles  or  by  prism  astigmatism  for  large  Lit¬ 
trow  angles.10  For  the  case  of  SHIMMER  in  which 
astigmatism  is  the  limiting  factor,  the  maximum  field 
of  view  can  be  obtained  by  use  of  a  prism  at  minimum 
deviation  with  angle  of  incidence  in  (see  Fig.  1)  given 
by 


( n 2  —  1)  2 n2  —  sin2  iq 

- 2 —  tan  in  — ^ —  =  tan  0L, 


sin2  p 


(4) 


where  n  is  the  prism  index  of  refraction.  Once  p  is 
determined  from  Eq.  (4),  the  prism  apex  angle  a  is 
given  by  the  condition  for  minimum  deviation 
through  the  prism: 


n  sin 


sin  in. 


(5) 


The  field-widening  prisms  also  introduce  resolution 
in  addition  to  the  grating  dispersion.  The  resulting 
resolving  power  R  is  given  by 


R=R0  1 


An  n  \ 
AX  n2  —  1/  ’ 


(6) 


where  the  correction  to  R0  is  due  to  dispersion  by  the 
prisms  with  index  of  refraction  n  at  wavelength  X. 


C.  Spatial  Imaging  Considerations 

Spatial  heterodyne  spectrometer  instruments  can  ob¬ 
tain  zero,  one,  or  two  dimensions  of  spatial  informa¬ 
tion  with  appropriate  choice  of  the  optics  prior  to  the 
interferometer.  When  the  scene  is  imaged  by  a  tele¬ 
scope  on  the  input  aperture,  each  point  on  the  inter¬ 
ferogram  collects  light  from  all  points  in  the  scene. 
In  this  case  the  interferogram  represents  the  average 
spectrum  over  the  field  of  view  of  the  telescope. 
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Fig.  2.  MAHRSI  detection  of  OH  at  70-km  tangent  height. 


However,  when  the  scene  is  imaged  on  the  localiza¬ 
tion  plane  of  the  fringes  near  the  gratings,  spatial  as 
well  as  spectral  information  can  be  obtained.  If  the 
scene  contains  only  one  dimension  of  spatial  informa¬ 
tion  as  is  the  case  with  SHIMMER  (we  expect  the  OH 
profiles  to  vary  spatially  only  perpendicular  to  the 
Earth’s  limb),  an  interferogram  for  each  altitude  can 
be  obtained  when  the  detector  is  read  out  row  by  row. 
In  situations  in  which  intensity  variations  parallel  to 
the  limb  distort  the  interferogram,  the  spatial  infor¬ 
mation  in  this  dimension  can  be  washed  out  by  use  of 
an  anamorphic  telescope.  We  can  obtain  two- 
dimensional  spatial  information  along  with  spectra 
by  scanning  the  field  of  view  during  an  observation 
while  time  tagging  detected  photons.11-12 

3.  SHIMMER  Instrument 

A.  Measurement  Considerations 

The  SHIMMER  instrument  has  been  specifically  de¬ 
signed  to  measure  OH  solar  resonance  fluorescence 
while  viewing  the  limb  from  low  Earth  orbit.  This 
measurement  is  a  rigorous  test  of  any  hyperspectral 
imaging  technology  because  retrieval  of  the  OH  radi¬ 
ance  profile  requires  highly  accurate  knowledge  of 
the  instrument’s  performance.  In  particular,  the  re¬ 
trieval  of  OH  spectra  from  limb  observations  entails 
identification  and  removal  of  the  bright  and  spec¬ 
trally  complex  ozone-attenuated  Rayleigh-scattered 
solar  background.  This  is  accomplished  when  we 
perform  a  nonlinear  least-squares  curve-fitting  pro¬ 
cedure  that  relies  on  high-spectral-resolution  mea¬ 
surements  to  properly  identify  and  separate  the 
spectral  components  of  the  limb  observations.  Di¬ 
rect  measurements  of  the  Moon  with  a  diffuser  are 
used  to  determine  the  solar  background  spectral 
shape  measured  by  SHIMMER  independent  of  ter¬ 
restrial  OH  emission. 

The  problem  is  illustrated  in  Fig.  2,  where  the  top 
panel  shows  a  MAHRSI-measured  OH  emission  spec¬ 
trum  at  70-km  tangent  height  and  the  underlying 
Rayleigh-scattered  solar  spectrum.  This  is  the  tan¬ 
gent  altitude  for  which  OH  is  most  prominent  against 


the  Rayleigh-scattered  solar  background.  The  lower 
panel  shows  the  OH  signal  extracted  when  the  Ray¬ 
leigh  background  is  subtracted  from  the  raw  spec¬ 
trum.  Because,  in  effect,  the  curve-fitting  algorithm 
subtracts  one  large  number  from  another  at  every 
wavelength,  leaving  the  residual  OH  spectrum,  any 
uncharacterized  and  uncorrected  instrumental  con¬ 
tributions  to  the  spectral  shape  will  result  in  large 
errors  in  the  retrieved  OH  intensity,  particularly  for 
observations  of  lower  tangent  heights.  The  error  in 
the  fit  of  the  background  shape  is  greatly  reduced 
when  we  utilize  the  full  bandwidth  of  the  spectral 
region  of  the  OH  fluorescence.  It  is  the  combination 
of  high  spectral  resolution  and  extended  bandwidth 
that  is  essential  to  the  accurate  identification  of  the 
background  spectrum  contribution.  Furthermore, 
as  the  altitude  intensity  profiles  are  inverted  to  yield 
OH  concentration  profiles,  the  intensity  errors  prop¬ 
agate  through  the  inversion  to  produce  even  larger 
concentration  uncertainties.  For  example,  large  re¬ 
trieval  and  inversion  errors  will  be  generated  if 
knowledge  of  the  instrument  function  and  its  varia¬ 
tion  with  wavelength  or  position  in  the  observed 
scene  is  incomplete  or  inaccurate,  if  the  spectral  re- 
sponsivity  of  the  instrument  is  not  well  character¬ 
ized,  or  if  off-axis  or  internally  scattered  light  is  not 
controlled  adequately.  The  design  specifications  for 
the  SHIMMER  instrument  are  summarized  in  Table 
1.  The  spectral  resolution,  passband,  field  of  view, 
and  responsivity  requirements  of  SHIMMER  were 
determined  by  the  previous  successful  OH  measure¬ 
ments  by  the  MAHRSI  instrument. 

The  optical  layout  is  shown  schematically  in  Fig.  3. 
The  500-mm  focal-length  refractive  telescope  focuses 
the  sky  on  the  gratings.  The  telescope  objective  di¬ 
ameter  limits  the  range  of  angles  incident  on  the 
gratings  while  also  imaging  the  scene  at  infinity  on 
the  grating  localization  plane.  This  property  will  be 
used  in  SHIMMER  to  image  altitude  along  the  col¬ 
umns  of  the  CCD  while  producing  one-dimensional 
interferograms  along  the  rows  while  limb  viewing 
from  space.  Telescope  aberrations  limit  the  altitude 
resolution  on  the  limb  to  0.87  km.  In  practice,  how¬ 
ever,  rows  will  be  coadded  on  the  chip  to  give  a  2.2-km 
(4.3-arc  min)  altitude  resolution. 

Out-of-band  rejection  is  achieved  by  a  6.6-cm- 
diameter  interference  filter  between  the  telescope 
and  the  interferometer  and  narrow-band  multilayer 
dielectric  coatings  on  the  three  fold  mirrors  in  the 
system.  All  transmitting  elements  have  been  anti¬ 
reflection  coated  for  308  nm.  The  interference  filter 
has  a  2.3-nm  FWHM  passband  centered  at  308.9  nm. 
Because  the  interferometer  in  this  design  cannot  dis¬ 
tinguish  between  wavelengths  at  +  AX  and  —  AX  about 
the  Littrow  wavelength,  the  violet  edge  of  the  filter 
passband  is  positioned  near  the  Littrow  wavelength 
(307.0  nm)  so  that  shorter  wavelengths  are  blocked. 
The  interferometer  uses  an  antireflection-coated, 
nonpolarizing  hexagonal  beam  splitter,  a  pair  of 
13.02°  apex  angle  field-widening  prisms,  and  a  pair  of 
1200  lines/mm,  75%  efficiency  gratings.  For  fused- 
silica  prisms,  1200-line/mm  gratings,  and  a  design 
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Table  1.  SHIMMER  Design  Specifications 


SHIMMER  PGSE 


Aspect  Attribute  Design 


System  Interference  filter 

2.3-nm  FWHM  centered 
at  308.9  nm 
307.0-310.5-nm  full 

width 

Entrance  optics 

500-mm  focal-length 

telescope 

Exit  optics 

7-element  relay  system 

Detector 

1034  X  1024  24  pm  X 

24  pan  CCD 

Gratings 

Clear  aperture 

20  mm  X  20  mm 

Groove  density 

1200  lines/mm 

Littrow  angle 

10.7° 

Littrow  wavelength 

307.0  nm 

Field-widening  prisms 

Clear  aperture 

22  mm  X  22  mm 

Wedge  angle 

13.02° 

Incident  angle 

8.73° 

Exit  angle 

10.7° 

Beam  splitter 

Clear  aperture 

28  mm  X  28  mm 

Field  of  view 

At  gratings 

10° 

On  sky 

2.3  X  2.3° 

Performance  Resolving  power 

53,500 

Spectral  resolution 

0.0058  nm 

Nonaliased  spectral 

2.95  nm 

range 

Achieved  spectral  range 

Filter  limited 

Sky  imaging 

1.7  arc  min 

Resources  Mass 

22  kg 

Volume 

47  cm  X  36  cm  X  24  cm 

Power 

27  W 

CCD  readout 

72  bytes/8  s 

Image  rate 

375  images/orbit 

Data  rate 

27  Mbytes/orbit 

wavelength  of  308.5  nm,  the  Littrow  angle  of  the 
grating  and  the  optimal  angle  of  refraction  of  the 
optical  axis  leaving  the  prism  are  nearly  equal  [see 
Eqs.  ( 1)— (3)  in  Section  2.  A];  hence  the  back  face  of  the 
prisms  and  the  gratings  are  nearly  parallel.  For 
ease  of  alignment  the  back  prism  face  and  the  grating 
face  were  made  parallel  for  SHIMMER  (see  Table  1). 
A  doublet  field  lens  and  a  multielement  relay  lens 
focus  the  plane  of  the  gratings,  which  is  the  fringe 
localization  plane,  on  the  CCD.  The  theoretical 
modulation  transfer  function  (MTF)  of  the  SHIM¬ 
MER  imaging  lenses  is  shown  in  Fig.  4. 

The  camera  used  in  SHIMMER  utilizes  a  thinned 
backilluminated  UV  antireflection-coated  CCD  man¬ 
ufactured  by  Scientific  Imaging  Technologies  Inc. 
with  approximately  58%  quantum  efficiency  at  308 
nm.  The  CCD  format  is  1024  X  1024  with  24  p,m  X 
24  |xm  pixels.  The  camera  also  includes  a  shutter 
mechanism,  thermoelectric  coolers,  cooling  fins,  and 
fans  to  convectively  dispel  heat  into  the  cabin  air. 

The  mechanical  design  for  SHIMMER  was  driven 
by  five  requirements:  (1)  locate  and  support  to  in¬ 


Fig.  3.  SHIMMER  block  diagram  and  optical  layout.  OCA,  or- 
biter  communication  adapter;  ISA,  industry-standard 
architecture-bus  computer;  PGSE,  payload  and  general  support 
electronics;  TEC,  thermoelectric  cooler. 


terferometric  precision  the  optical  elements  of  the 
interferometer  in  a  manner  suitable  for  Space  Shut¬ 
tle  launch  and  landing  loads;  (2)  locate  and  fasten  the 
optical  components  shown  in  Fig.  3  to  ensure  accu¬ 
rate  alignment;  (3)  constrain  the  mechanical  enve¬ 
lope  of  the  instrument  to  the  interior  dimensions  of 
the  Space  Shuttle  middeck  lockers;  (4)  provide  a  sim¬ 
ple,  robust  mating  assembly  to  locate  and  fasten  the 
instrument  to  the  Space  Shuttle  crew  hatch;  and  (5) 
meet  all  manned-flight  program  safety  requirements. 

The  interferometric  elements,  i.e.,  the  beam  split¬ 
ter,  prisms,  and  gratings,  are  mounted  in  a  VascoMax 
steel  structure  shown  in  Fig.  5  and  labeled  inter¬ 
ferometer  in  Fig.  6.  The  beam  splitter  is  positioned 
below  the  triangular  cap  visible  at  the  top,  and  the 
input  and  output  windows  are  on  the  backside.  The 


Fig.  4.  Theoretical  exit  optics  MTF.  The  three  solid  curves  are 
for  on-axis,  70%  field  and  full  field  points  on  the  gratings.  The 
dashed  curve  shows  the  diffraction-limited  MTF.  The  pixels  on 
the  CCD  limit  the  maximum  frequency  to  21  cycles/mm. 
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Fig.  5.  SHIMMER  components  on  the  flight  baseplate  in  the 
laboratory. 


design  allows  precise  and  stable  positioning  of  the 
prisms  and  gratings  relative  to  the  beam  splitter. 
The  prism  holders  are  attached  to  the  grating  hold¬ 
ers,  and  a  spring-loaded  turret  assembly  holds  the 
grating  holders  tightly  to  a  three-point  kinematic 
mount  integral  to  the  cell  structure.  An  interior 
chamber  encloses  the  CCD  camera  and  is  walled  off 
by  light-tight  bulkheads  and  baffle  tubes.  Both  air¬ 
flow  apertures  open  into  this  isolated  space.  Each 
aperture  is  covered  by  a  50-p.m  mesh  screen.  The 
entrance  aperture  is  protected  by  a  transparent 
Lexan  cover  that  will  be  removed  on  orbit  by  the 
flight  crew.  The  cover  plate  of  the  housing  unit  also 
provides  attach  points  for  lifting  the  handles. 

B.  Operational  Considerations 

As  described  above,  SHIMMER  was  designed  to  con¬ 
duct  a  spaceflight  proof  test  of  the  application  of  the 
SHS  technique  to  remote  sensing  of  the  middle  at¬ 
mosphere.  The  test  will  be  to  observe  the  vertical 
profile  of  OH  solar  resonance  fluorescence  when  the 
Earth’s  Sun-illuminated  horizon  is  imaged  through 
the  UV-transmitting  optical  window  in  the  Space 
Shuttle  orbiter  crew  hatch.  Tests  on  three  Shuttle 
windows  have  shown  that  the  panes  are  not  flat 


Wilson  Stud  Assemblies 


enough  to  cause  additional  interference  effects  in 
SHIMMER  interferograms.  The  instrument  is  at¬ 
tached  to  the  orbiter  crew  hatch  with  a  bracket  as¬ 
sembly  that  utilizes  a  circle  of  threaded  studs  located 
around  the  optical  window.  As  shown  in  Fig.  7,  the 
bracket  includes  a  hinged  channel  that  supports  the 
instrument  and  a  power  converter  module.  The  as¬ 
sembly  is  hinged  to  permit  folded  stowage  in  a  mid¬ 
deck  locker. 

The  crew  will  maneuver  the  orbiter  to  an  attitude 
that  will  place  the  SHIMMER  line-of-sight  tangent  to 
the  limb  at  a  nominal  altitude  of  65  km  while  placing 
the  window  out  of  the  Sun’s  direct  illumination.  The 
projected  2.3  deg  instrument  field  of  view  will  then 
have  an  upper  boundary  at  100-km  altitude  and  a 
lower  boundary  at  30-km  altitude.  At  a  wavelength 
of  308  nm,  the  vertical  radiance  profile  of  Rayleigh- 
scattered  sunlight  has  a  well-defined  maximum  at 
approximately  41-km  altitude.  The  flight  crew  will 
confirm  the  pointing  of  the  instrument  by  observing 
the  placement  of  this  radiance  peak  in  the  field  of 
view  displayed  on  a  laptop  computer.  The  crew  then 
will  initiate  several  sequences  of  nominal  (1.2-s)  ex¬ 
posures.  The  data  will  be  downlinked  in  near  real 
time  to  an  operations  control  center  at  Johnson  Space 
Center  where  the  SHIMMER  science  team  will  ana¬ 
lyze  the  data  and  provide  feedback  to  the  crew.  Af¬ 
ter  the  flight,  the  data  will  be  analyzed  to  produce  OH 
number  density  profiles  by  use  of  algorithms  vali¬ 
dated  in  the  laboratory  and  on  MAHRSI  data. 

4.  Instrument  Characterization  in  the  Laboratory 

An  extensive  series  of  laboratory  tests  have  verified 
that  SHIMMER  has  met  design  expectations.  The 
most  basic  test  is  to  observe  the  interferogram  pro¬ 
duced  by  a  monochromatic  source  within  the  instru¬ 
ment’s  spectral  filter  passband  ( %T  >  35%  from  307.8 
to  309.6  nm).  A  Zn  pen-ray  lamp  with  a  single 
bright  emission  at  307.6  nm  was  used  for  optimiza¬ 
tion  and  testing.  The  spectral  width  of  the  emission 
is  not  resolved  by  the  interferometer.  The  fringe 
pattern  recorded  by  SHIMMER  with  the  Zn  pen-ray 
source  is  shown  in  Fig.  8.  Note  that  all  spectral 
information  is  in  the  horizontal  dimension,  and  the 
vertical  dimension  is  exclusively  an  image  of  the  ver¬ 
tical  intensity  distribution  of  the  source.  When  ap- 
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Fig.  8.  Interferogram  recorded  in  the  laboratory. 


plied  to  the  OH  limb-viewing  problem,  the  vertical 
dimension  will  be  aligned  perpendicular  to  the  limb 
so  that  each  row  of  the  CCD  will  record  the  spectrum 
observed  at  distinct  viewing  angles  corresponding  to 
distinct  tangent  heights. 

As  discussed  in  Section  2.A,  grating  figure  errors 
and  interferometer  refractive-index  inhomogeneities 
distort  the  Fizeau  fringe  pattern  recorded  at  the 
CCD,  but  with  the  high-quality  optics  of  SHIMMER, 
these  distortions  are  negligible.  The  interferogram 
was  corrected  for  dark  field  and  nonuniformity,  then 
Fourier  transformed  to  produce  the  spectrum  shown 
in  Fig.  9.  Because  the  Zn  line  at  307.6  nm  is  unre¬ 
solved  by  the  instrument,  this  spectrum  represents 
the  achieved  instrument  function. 

The  passband  of  the  instrument  is  determined  by 
the  spectral  filter,  the  spectral  resolution,  the  number 
of  pixels  sampling  the  interferogram,  and  the  image 
quality  at  the  CCD.  In  these  tests,  we  focused  the 
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Fig.  9.  Measured  Zn  emission  feature  at  307.6  nm. 


image  of  the  gratings  onto  1024  pixels,  so  the  number 
of  spectral  elements  available  without  aliasing  is  512. 
Assuming  perfect  imaging,  the  maximum  passband 
should  be  512  times  the  spectral  resolution  (0.0058 
nm  as  verified  by  means  of  OH  emission  as  described 
below),  or  approximately  3.0  nm. 

Tests  with  a  polychromatic  source  were  conducted 
with  a  microwave-excited  OH  lamp.  These  tests 
verified  that  the  spectral  resolution,  bandwidth,  and 
field  of  view  all  meet  design  specifications.  With  the 
gratings  adjusted  for  a  Littrow  wavelength  of  306.42 
nm,  the  OH  spectrum  shown  in  the  first  panel  of  Fig. 
10  was  recorded.  We  generated  the  spectrum  in  the 
second  panel  using  OH  g  factors  assuming  a  solar 
source  and  a  rotational  temperature  of  2000  K.  Fig¬ 
ure  10  shows  that  the  two  features  at  307.526  and 
308.607  nm  are  recorded  in  spectral  resolution  ele¬ 
ments  213  and  395.  This  coupled  with  the  one-pixel¬ 
wide  line  shape  in  the  Zn  spectrum  in  Fig.  9  verify  a 
spectral  resolution  of  0.0059  nm,  or  a  resolving  power 
of  53,500  as  designed.  The  second  panel  shows  the 
OH  spectrum  and  filter  transmission  folded  about  the 
Littrow  wavelength  of  306.42  nm.  As  mentioned 
above,  the  interferometer  cannot  distinguish  between 
wavelengths  at  +AX  and  —AX  about  the  Littrow 
wavelength,  but,  because  the  short-wavelength  edge 
of  the  filter  passband  is  near  the  Littrow  wavelength, 
the  wavelengths  to  the  violet  of  Littrow  are  blocked. 

We  measured  the  average  spectral  responsivity  of 
the  instrument  using  a  calibrated  Labsphere  Inc.  in¬ 
tegrating  sphere  with  an  8-in.-  (20-cm-)  diameter  ap¬ 
erture,  internal  Spectralon  diffuse  reflecting  tiles, 
and  multiple  halogen  lamp  sources.  The  absolute 
spectral  radiance  of  the  sphere  was  measured  at  the 
NASA  Goddard  Space  Flight  Center  Laboratory  for 
Atmospheres  in  collaboration  with  the  Shuttle  Solar 
Backscatter  Ultraviolet  experiment  team.13-14  The 
SHIMMER  field-of-view  footprint  is  within  the 
sphere’s  aperture.  The  total  signal  at  the  CCD  (dig¬ 
ital  number,  DN)  is  related  to  the  total  inband  photon 
flux  at  the  SHIMMER  input  aperture  by 

S(DN)  =  4>(phot/s)Tint(s)T'filtT,optQECCD(e  /phot) 

x  G(DN/0,  (7) 

where  Tint  is  the  exposure  time,  Tmt  and  T'opt  are  the 
average  spectral  filter  and  optical  element  transmis¬ 
sion  over  the  passband  of  the  instrument,  and  QECCD 
and  G  are  the  quantum  efficiency  and  gain  of  the 
CCD.  The  total  inband  photon  flux  cf>  is  related  to 
the  spectral  radiance  of  the  integrating  sphere  Lx  [in 
kiloRayleighs  (kR)  per  nanometer]  by 

109 

4>(phot/s)  =  Lx  (kR/nm)  - —  AX(nm)A(cm2)H(sr), 


where  AX  is  the  bandwidth,  A  is  the  input  aperture 
area,  and  O  is  the  field  of  view.  The  measured  re¬ 
sponsivity  of  SHIMMER  is  R  =  4.34  X  104  DN  s~4 
kR-1,  the  derived  optical  transmission  is  Topt  =  0.18, 


1  March  2002  /  Vol.  41,  No.  7  /  APPLIED  OPTICS 


1349 


SHIMMER  OH  Measurement 


ui 

c 

CD 

E 

c 

cm 

in 

o 

o 

o 


CL 


1x10 
8x10 
6x1  0 
4x10 
2  x  1  0 


100 


2G0 


300 


400 


306.5  307.0  307.5  308.0  308.5  309.0 

Wavelength  (nm) 

Fig.  10.  OH  spectrum  measured  in  the  laboratory  compared  with  the  theoretical  OH  spectrum. 


and  the  calculated  throughput  is  ToptTfiltAfl  =  4.5  X 
10 ~3,  all  meeting  instrument  design  goals. 

5.  SHIMMER  Earth  Limb  Measurement  Simulation 

Model  interferograms  were  generated  to  predict  the 
performance  of  the  instrument  while  viewing  the 
limb  from  space.  High-resolution  OH  g  factors,  so¬ 
lar  spectrum,  and  spectral  filter  contributions  are 
combined  to  produce  an  effective  incident  spectrum. 
The  CCD  rows  representing  a  1.2-s  integration  were 
binned  into  32  row  slices  to  increase  the  signal- 
to-noise  ratio,  with  the  OH  and  solar  contributions 
appropriately  weighted  based  on  MAHRSI  measure¬ 
ments  of  OH  and  Rayleigh-scattered  background  in¬ 
tensity  collected  during  the  STS-85  Shuttle  mission. 
Assuming  a  unit  power  image  of  the  gratings  at  the 
CCD,  the  DNs  per  pixel  per  slice  [DN(pixel/slice)] 
that  are  due  to  altitude-dependent  spectral  radiances 
L0 h  and  Lsun  (kR)  are  calculated  with  the  equation, 

(VO'+l)  fzti+l)  C" 

‘S’ccdIC 7)  =  Lou(cr,  y)R(cr) 

Jy(j)  Jx(i)  Jo 

X  (1  +  cos{2tt[4(ct  —  cr q) x  tan  0]})do-dxdy 

(M./+1)  f” 

+  Lsun(a,  y)fi(a) 

Jy(j)  Jx(i)  Jo 

X  (1  +  cos{2tt[4(ct  —  cr0)x  tan  0]})do-dxdy 

DN(pixel/slice),  (9) 

where  i  is  the  CCD  pixel  number  along  rows;./  is  the 
slice  number  corresponding  to  one  of  32  tangent 
height  altitudes  in  the  range  30-94  km;  x(i)  and  y(j) 
are  the  location  on  the  grating  imaged  onto  the  corner 
of  the  CCD  region  (i,j);  and  the  responsivity  R((t)  that 
incorporates  entrance  aperture  size,  pixel  field  of 
view,  optical  transmission,  CCD  quantum  efficiency, 


integration  time,  and  analog-to-digital  converter  gain 
converts  radiance  at  the  limb  (kR)  to  signal  readout 
from  the  CCD  camera  [DN(pixel/slice)].  Figure  11 
shows  the  interferogram  produced  by  anticipated  OH 
emissions  and  Rayleigh-scattered  sunlight.  The 
emissions  from  lower  altitudes,  approximately  a  fac¬ 
tor  of  40  brighter  than  those  from  higher  altitudes, 
are  imaged  at  the  bottom  of  the  CCD.  The  line  at 
zero  optical  path  difference,  at  which  the  frequency 
components  that  are  due  to  all  wavelengths  and  in¬ 
cidence  angles  are  in  phase,  is  evident  in  the  central 
column.  The  upper  panel  shows  the  signal  (solid 
curves)  and  noise  (dashed  curves)  components 
present  in  the  interferogram  at  the  CCD  that  are  due 
to  atmospheric  emissions  and  instrumental  effects. 
All  the  curves  shown  in  the  upper  panel  are  in  units 
of  DN(pixel/slice).  These  curves  indicate  the  rela¬ 
tive  contribution  to  the  interferogram  as  a  function  of 
altitude  for  each  of  the  signal  and  noise  sources. 
Note  that  the  total  signal  and  noise  are  dominated  by 
Rayleigh-scattered  background  at  all  altitudes. 
Note  also  that  instrumental  noise  contributions  are 
only  significant  above  approximately  85  km  and  that 
the  total  noise  becomes  comparable  to  or  dominates 
the  OH  signal  only  above  85  km  and  below  40  km. 

We  made  estimates  of  the  signal-to-noise  ratio  at 
each  altitude  by  retrieving  OH  intensity  as  a  function 
of  altitude  from  the  interferogram  in  the  same  way  as 
will  be  done  during  the  mission.  In  this  procedure 
we  subtract  the  dark-field,  flat-field  calibration  to  ra¬ 
diance  units;  Fourier  transform,  phase  correct  and 
perform  a  nonlinear  curve  fitting  to  remove  the 
Rayleigh-scattered  background;  and  integrate  over 
the  OH  features.  The  bottom  panel  shows  the  re¬ 
trieved  OH  spectral  features  (horizontal)  as  a  func¬ 
tion  of  altitude  (vertical).  The  retrieved  radiance  as 
a  function  of  altitude  (shown  in  Fig.  12)  is  the  same  as 
the  radiance  profile  input  to  the  instrument  model  in 
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Fig.  11.  (upper  panel)  Modeled  signal  and  noise  contributions,  (middle  panel)  1024  X  32  pixel  interferogram,  and  (lower  panel)  512  X  32 
spectral  resolution  element  retrieved  OH. 


the  first  place,  as  it  should  be.  The  estimated  pre¬ 
cision  of  the  measurement  as  a  function  of  altitude 
(shown  in  Fig.  13)  is  determined  not  only  by  the 
signal-to-noise  ratio  in  the  interferogram,  but  also  by 
the  bandwidth  and  spectral  shape  of  the  incident 
spectrum. 

6.  Future  Developments 

We  are  currently  pursuing  further  miniaturization  of 
an  OH  spatial  heterodyne  spectrometer  with  the  sup¬ 
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Fig.  12.  Modeled  OH  altitude  intensity  profile. 


port  of  NASA’s  Planetary  Instrument  Definition  and 
Development  Program.  Further  reduction  in  in¬ 
strument  size  and  weight  will  be  achieved  when  the 
interferometer  components  are  optically  contacted  to 
form  a  monolithic  interferometer.  This  technique 
makes  the  complex  VascoMax  structure,  which  is 
used  in  SHIMMER  to  hold  the  optical  components  in 
precise  proximity,  obsolete.  Furthermore,  the  align¬ 
ment  of  a  monolithic  interferometer  is  virtually  im¬ 
pervious  to  vibrations,  resulting  in  a  major  risk 
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Fig.  13.  Estimated  OH  radiance  measurement  precision. 
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Fig.  14.  Full-scale  model  of  a  miniature  monolithic  version  of 
SHIMMER. 


reduction  for  spaceflight  applications.  Figure  14 
shows  a  photograph  of  the  monolithic  instrument 
presently  under  construction.  It  will  have  200  times 
the  throughput  of  MAHRSI  and  a  spectral  resolution 
of  0.012  nm.  We  plan  to  report  on  the  performance 
of  the  monolithic  instrument  in  the  near  future. 

The  latest  achievements  in  SHS  and  the  ongoing 
research  confirm  that  this  technique  is  remarkably 
well  suited  for  the  investigation  of  planetary  atmo¬ 
spheres  from  space.  For  example,  longstanding  is¬ 
sues  such  as  the  polar  mesospheric  cloud 
phenomenon  on  Earth  or  the  detailed  understanding 
of  the  photochemistry  in  the  Martian  atmosphere  can 
be  addressed  with  global  [OH],  [NO],  and  tempera¬ 
tures  provided  by  high-spectral-resolution  UV  limb 
spectra  acquired  from  an  orbiting  platform. 
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